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Novel relativistic expressions are used to calculate the weak decay constants of 
pseudoscalar and vector mesons within the constituent quark model. Meson wave 
functions satisfy the quasipotential equation with the complete relativistic potential, 
pgj ■ New contributions, coming from the negative-energy quark states, are substantial 

for the light mesons, significantly decrease the values of their decay constants and, 
thus, bring them into agreement with experiment. For heavy-light mesons these 
contribution are much less pronounced, but permit to reduce uncertainties of the 
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Q_i' predicted decay constants. Their values agree with the results of lattice calculations 

and experimental data. 

<D : 

' PACS numbers: 13.20.-v, 14.40.-n, 12.39.Ki 

• • i 

> ■ 

The weak decay constants of pseudoscalar and vector mesons belong to their most impor- 
tant characteristics, which enter in various decay rates. Many efforts were undertaken to cal- 
culate these constants within lattice QCD (both quenched and unquenched) E l S P, HI IH El , 
QCD sum rules 0, 0, EJ, and constituent quark models 0, Q E El Q El Q- At 
present, the decay constants of light mesons are measured with high precision, while in the 
heavy-light meson sector only D and D s meson decay constants are available with rather 
large errors [l7| . Recently, a relatively precise experimental value for the D meson decay 
constant was presented by the CLEO Collaboration [18|. Therefore it is actual to recon- 
sider the meson decay constants treating quarks, composing the meson, in a consistently 
relativistic way. Such procedure was formulated and successfully applied for light mesons in 
the papers (19]. In this letter we evaluate new contributions to relativistic expressions for 
the meson decay constants coming from the negative-energy quark states both for light and 
heavy-light mesons. We use the meson wave functions satisfying the quasipotential equation 
with the complete relativistic potential in order to obtain new, more accurate predictions 
for the meson decay constants. 

In the quasipotential approach a meson is described by the wave function of the bound 
quark- ant iquark state [20], which satisfies the quasipotential equation of the Schrodinger 
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type 

(b 2 (M) p 2 \ /■ c/ 3 g 



\ 2^7 i (2tt 

where the relativistic reduced mass is 



*m(p) = / 7^x^(P> q; M)* v (q), (i) 



£^2 M 4 - (m 2 - ml) 2 



and Ei, E 2 are given by 



M 2 — m 2 , + m? M 2 — m? + m 2 , 
E 1 = 2 — — 1 -, E 2 = ±— — 2 -. 3 

Here M = E\ + E 2 is the meson mass, mi^ are the quark masses, and p is their relative 
momentum. In the center-of-mass system the relative momentum squared on mass shell 
reads 

2 = [M 2 -( mi + m 2 ) 2 ][M 2 -( mi -m 2 ) 2 ] 
1 ' AM 2 ' 1 1 

The kernel V(p, q; M) in Eq. is the quasipotential operator of the quark- ant iquark in- 
teraction. It is constructed with the help of the off-mass-shell scattering amplitude, projected 
onto the positive energy states. Constructing the quasipotential of the quark-antiquark in- 
teraction, we have assumed that the effective interaction is the sum of the usual one-gluon 
exchange term with the mixture of long-range vector and scalar linear confining potentials, 
where the vector confining potential contains the Pauli interaction. The quasipotential is 
then defined by 1 

V(p, q; M) = ui(p)u 2 (-p)V(p, q; M)ui(q)u a (-q), (5) 

with 

V(p, q; M) = V(p - q) = ^a s D, u (k)^Y 2 + (k^I^ + (k), 

where a s is the QCD coupling constant, D^ v is the gluon propagator in the Coulomb gauge 
and k = p — q; 7 M and u(p) are Dirac matrices and spinors. The effective long-range vector 
vertex is given by 

rv(k) = 7m + 7^<V r > ( 6 ) 

where k is the Pauli interaction constant characterizing the anomalous chromomagnetic 
moment of quarks. Vector and scalar confining potentials in the nonrelativistic limit reduce 
to 

VcoAr) = (l-s)(Ar + B), 

Konf(r-) = e(Ar + B), (7) 

reproducing 

KoBf(r) = (r) + VLt(r) = Ar + B, (8) 



1 In our notation, where strong annihilation processes are neglected, antiparticles are described by usual 
spinors taking into account the proper quark charges. 
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TABLE I: Masses of the ground state light and heavy-light 


mesons (in MeV). 


Meson 


^~theor 


M cxp pDG ^ 


7T 


154 


139.57 


P 


776 


775.8(5) 


K 


482 


493.677(16) 


K* 


897 


891.66(26) 


<t> 


1038 


1019.46(2) 


D 


1872 


1869.4(5) 


D* 


2009 


2010.0(5) 


D s 


1967 


1968.3(5) 


D* s 


2112 


2112.1(7) 


B 


5275 


5279.0(5) 


B* 


5326 


5325.0(6) 


B s 


5362 


5369.6(2.4) 


IK 


5414 


5416.6(3.5) 



where e is the mixing coefficient. 

All the model parameters have the same values as in our previous papers [13, EJ. The 
constituent quark masses m u = = 0.33 GeV, m s = 0.5 GeV, m c = 1.55 GeV, m& = 4.88 
GeV and the parameters of the linear potential A = 0.18 GeV 2 and B = —0.3 GeV have 
the usual values of quark models. The values of the mixing coefficient of vector and scalar 
confining potentials e — — 1 and the universal Pauli interaction constant k — — 1 are specific 
for our model. 

The quasipotential (0) can be used for arbitrary quark masses. The substitution of the 
Dirac spinors into (0) results in an extremely nonlocal potential in the configuration space. 
Clearly, it is very hard to deal with such potentials without any additional transformations. 
In oder to simplify the relativistic qq potential, we make the following replacement in the 
Dirac spinors [19j, |2l| : 

ei, 2 (p) = + P 2 -> E 1>2 . (9) 

This substitution makes the Fourier transformation of the potential (0) local, but the re- 
sulting relativistic potential becomes dependent on the meson mass in a very complicated 
nonlinear way. We consider only the meson ground states, which further simplifies our 
analysis, since all terms containing orbital momentum vanish. The detailed expressions for 
the relativistic quark potential can be found in Ref. [19J. Here we use these formulas for 
the calculation of the ground state meson masses. We solve numerically the quasipotential 
equation with the local fully relativistic potential, which includes both spin-independent and 
spin-dependent parts. As a result we get the relativistic wave functions of the ground state 
mesons which depend nonperturbatively on the meson spin (i.e. the pseudoscalar and vector 
meson wave functions are different). These wave functions are used below for calculating 
the decay constants of light and heavy mesons. The obtained masses of the pseudoscalar 
and vector mesons are given in Table in comparison with the experimental data ^?|- The 
overall good agreement of our predictions with experiment is found. 

The decay constants fp and fy of the pseudoscalar (P) and vector (V) mesons parame- 
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*(M,p) 





(1) 



(2) 



(3) 



(4) 



FIG. 1: Weak annihilation diagram of the meson. Solid and bold lines denote the positive- and 
negative-energy part of the quark propagator, respectively. Dashed lines represent the interaction 
operator V. Dashed ovals depict the projected wave function ^mk(p)- 



terize the matrix elements of the weak current Jjf = q\J^q2 = <?i 7 ^(l — 75)^2 between the 
corresponding meson and vacuum states. They are denned by 



(0|gi7 M 7 5 g 2 |P(K)) 

(0|g l7 W(K,£)> 



fvMve", 



(10) 
(11) 



where K is the meson momentum, e M and My are the polarization vector and mass of the 
vector meson. This matrix element can be expressed through the two-particle Bethe-Salpeter 
wave function \1/(M, p) in the quark loop integral (see Fig. 



I ^^H(1-7 5 )*(M,P)} 



(27T) 



(12) 



where the trace is taken over spin indices. Integration over p° in Eq. (|T2)l allows one to pass 
to the Fourier transform of the single-time wave function in the meson rest frame 



tt(M,p) 



2^ 



*(M,p). 



(13) 



This wave function contains both positive- and negative-energy quark states. Since in the 
quasipotential approach we use the wave function 1 Pmk(p) projected onto the positive- 
energy states, it is necessary to include additional terms which account for the contributions 
of negative-energy intermediate states. Within perturbation theory the weak matrix element 
(112)1 is schematically presented in Fig.0 The first diagram in the right hand side corresponds 
to the simple replacing of the wave function (JT3)l ty(M, p) by the projected one \I/a/k(p)- 2 
The second and third diagrams account for negative-energy contributions to the first and 
second quark propagators, respectively. The last diagram corresponds to negative-energy 
contributions from both quark propagators. 

Thus in the quasipotential approach this decay amplitude has the form 



0\JW\M(K] 



-jT^y MPl) J™ «2(P 2 )*MK(P) + 

(-)/„/ ^,0 tW a (+) 



V2M 



d 3 pd 3 p' _ 



ArVih^A^ (P' 2 )7 C 
M + etip') -e 2 (p') 



r 2 M 2 (P2)^(p 



(2tt) 6 

p')*a/k(p) + (1 <-> 2) 



2 The contributions with the exchange by the effective interaction potential V which contain only positive- 
energy intermediate states are automatically accounted for by the wave function itself. 
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rd 3 pd 3 P '_ Ai- ) (pl)7 ^- ) (^) 7 V , , .1 

(14) 

where p[^ 2 = K/2 ± pO; e(p) = Vp 2 + m 2 ; matrices r 12 denote the Dirac structure of the 
interaction potential (JHJ) for the first and second quark, respectively, and thus rir 2 V(p — 
p') = V(p — p'). The factor y2M follows from the normalization of the quasipotential wave 
function. The positive- and negative-energy projectors have standard definition 

A(±)( p ) = ± (^7° + 7°(7P)) 
2e(p) 

The ground-state wave function in the rest frame of the decaying meson \I/ M (p) = v I / a/o(p) 
can be expressed through a product of radial $m(p), spin Xss' and colour 4> qiq2 wave functions 

*Af(p) = $ m (p);w</w ( 15 ) 

Now the decay constants can be presented in the following form 

f _ , f(2+3) , f (4) 

where the terms on the right hand side originate from the corresponding diagrams in Fig. ^ 
and parameterize respective terms in Eq. ()14J) . In the literature (loL 14j usually only the 



first term is taken into account, since it provides the nonrelativistic limit, while other terms 
give only relativistic corrections and thus vanish in this limit. Such approximation can be 
justified for mesons containing only heavy quarks. However, as it will be shown below, for 
mesons with light quarks, especially for light mesons, other terms become equally important 
and their account is crucial for getting the results in agreement with experimental data. 

The matrix element (|14j) and thus the decay constants can be calculated in an arbitrary 
frame and from any component of the weak current Such calculation can be most easily 
performed in the rest frame of the decaying meson from the zero component of the current. 
The same results will be obtained from the vector component; however, this calculation is 
more cumbersome, since then the rest frame cannot be used and, thus, it is necessary to 
take into account the relativistic transformation of the meson wave function from the rest 
frame to the moving one with the momentum K. It is also possible to perform calculations 
in an explicitly covariant way using methods proposed in 2^ ■ 



The resulting expressions for decay constants are given by 



\ 1/2 

• - , , , < i (/') • n>:\ ( e 2 (p) + in-, 

Jpy 




26i (p) 
P 2 



ei(p) +m 1 ][e 2 (p) + m 2 ] 




(17) 



f(2+3) 

Jpy 




'ej(p) 



(p)[ei(p) + rni 



V 2e 2 (p) 
etjp) +mi ' 

e 2 (p) +m 2 . 



m 2 \ 



1/2 



"M-ei(p) -€ 2 (p) 
_M + e 1 (p) -e 2 (p) 

-2) Wp), 



(18) 
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TABLE II: Different contributions to the pseudoscalar and vector decay constants of light and 



heavy mesons 


(in MeV). 


The notations 


are taken according 


to Eqs. (HIIJ) and lf2"U|). 




Constant 


f NR 
JM 


f (i) 

J M 


, (2+3) ,(4) 
JM ~i~JM 


,,(2+3) ,(4)w f (l) 

Um "t Hi )/ Jm 


fhi 


u 


1290 


515 


-391 


-76% 


124 


n 


490 


402 


-183 


-46% 


219 


fx 


783 


353 


-198 


-56% 


155 


fK* 


508 


410 


-174 


-42% 


236 


u 


511 


415 


-170 


-41% 


245 


fD 


376 


275 


-41 


-15% 


234 


fD* 


391 


334 


-24 


-7% 


310 


fD s 


436 


306 


-38 


-12% 


268 


fD* 


447 


367 


-52 


-14% 


315 


fB 


259 


210 


-21 


-10% 


189 


fB* 


280 


235 


-16 


-7% 


219 


fB s 


300 


238 


-20 


-8% 


218 


fB* s 


316 


264 


-13 


-5% 


251 



f (4) 

Jpy 



12 f d 3 p ( ei(p) + mi 
, 2ei(p) , 
P 2 




1/2 



1/2 

e 2 {p)+m 2 \ M - ei(p) - e 2 (p) 



2e 2 (p) 



M + ei(p) + e 2 (p) 



ei(p) +mi][e 2 (p) + m 2 ] 

2 



P 



e?(p)^(p) ( e ib) + m i][ e 2{p) + m 2 ] 



*p,y(p), 



(19) 



with Xp = —1 and Ay = 1/3. Here e is the mixing coefficient of scalar and vector confining 
potentials and the long-range anomalous chromomagnetic quark moment K (JB|) is put 
equal to —1. Note that fp vanishes for pseudoscalar mesons with equal quark masses, 
such as the pion. The positive-ener gy contribution (JT7j) reproduces the previously known 
expressions for the decay constants [lfj. The negative-energy contributions (JTSj) and (|T9"j) 
are new and play a significant role for light mesons. 

In the nonrelativistic limit p 2 /m 2 — > the expression (fTTj) for decay constants gives the 
well-known formula 



NR 

py 



12 



M, 



I* 



py 



py 



(20) 



0. All other contributions vanish 



where \l/py(0) is the meson wave function at the origin r 
in the nonrelativistic limit. 

In Table ILT1 we present our predictions for the light and heavy-light meson decay constants 
calculated using the meson wave functions which were obtained as the numerical solutions 
of the quasipotential equation. 3 The values 4 of /^ R , obtained from the nonrelativistic 



3 We roughly estimate the uncertainties in our calculations to be about ten MeV for light mesons and of a 
several MeV for heavy-light mesons. 

4 For the evaluation of the relativistic wave functions were used. Thus the difference of the pseudoscalar 
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TABLE III: Pseudoscalar and vector decay constants of light mesons (in MeV). 



Constant this work [10] [11] [12] [16] Lattice [2] Lattice [3] Experiment [17] 





124 


180 


131 


219 138 


126.6(6.4) 


129.5(3.6) 


130.70(10)(36) 




155 


232 


155 


238 160 


152.0(6.1) 


156.6(3.7) 


159.80(1.4) (44) 












239.4(7.3) 




f 220(2)* 




219 


220 


207 


238 




\ 209(4)** 




236 


267 


241 


241 


255.5(6.5) 




217(5)t 


fa 


245 


336 


259 




270.8(6.5) 




229(3)* 



* derived from the experimental value for T p o^ e + e -. 
** derived from the experimental value for Y T ^ pVr . 

< derived from the experimental value for T t ^,k*u t - 

* derived from the experimental value for r^ e + e - . 



expression (}2T)j) . as well as the values of different contributions /^' 2 ' 3 ' 4 ^ ffTTj) - (fE?j) and the 
complete relativistic values of fu (QUI) are given. In Table IIHI we compare our results 5 for 
the decay constants Jm of light mesons with other quark model predictions [l(J EH EH G3 , 
recent values from two- and three-flavour lattice QCD and available experimental data 
It is clearly seen that the nonrelativistic predictions are significantly overestimating 
all decay constants, especially for the pion (almost by a factor of 10). The partial account 
of relativistic corrections by keeping in Eq. (fTo|) only the first term fj^ (fTTj) . which is 
usually used for semirelativistic calculations, does not substantially improve the situation. 
The disagreement is still large. This is connected with the anomalously small masses of 
light pseudoscalar mesons exhibiting their chiral nature. In the semirelativistic quark model 
[l0( the pseudoscalar meson mass is replaced by the so-called mock mass Mp, which is 
equal to the mean total energy of free quarks in a meson, and with our wave functions: 
M n = 2(e,(p)> « 1070 MeV (~ 8M W ) and M K = (e q (p)) + (e s (p)) w 1232 MeV (~ 2.5M K ). 

Such replacement yields values of fp which are still ~ 1.4 times larger than experimental 
ones (cf. 0). As we see from Table |Hj it is not justified to neglect contributions of the 
negative energy intermediate states for light meson decay constants. Indeed, the values of 
/m +3 ^ + Im are large and negative (reaching —76% of for the pion) thus compensating 
the overestimation of decay constants by the positive-energy contribution fff. This is the 
consequence of the smallness of the light pseudoscalar meson masses compared to the energies 
of their constituents. The negative-energy contributions (|18|) . (|19j) are proportional to the 
ratio of the meson binding energy M — €i(p) — € 2 (p) to its mass and quark energies. For 
mesons with heavy quarks this factor is small and leads to the suppression of negative-energy 
contributions. This results in the dominance of the positive-energy term ■ Indeed the 
negative-energy terms for heavy-light D and B mesons give 10 — 15% contributions (see 



and vector decay constants in this limit results from the difference of the corresponding relativistic wave 
functions. 

5 In our model p and u mesons are degenerate, therefore their decay constants are equal. The experimental 
value for the decay constant of the ui meson, derived from r w ^ e + e - ^tJj is = 195(3) MeV. 



8 



TABLE IV: Pseudoscalar and vector decay constants of heavy mesons (in MeV). 
Constant Quark models Lattice QCD QCD sum rules Experiment 





this work 


[15] 


m 


[4,5] 


W rs] [9] 


fD 


234 


230(25) 


225(14)(40) 


201(3)(17) 


203(20) 195(20) 




268 


248(27) 


267(13)(48) 


249(3)(16) 


235(24) 


fD a 1 Id 


1.15 


1.08(1) 




1.24(1)(7) 


1.15(4) 


fa 


189 


196(29) 


208(10)(29) 


216(9)(19)(6) 203(23) 206(20) 210(19) 


Ib s 


218 


216(32) 


250(10)(35) 


259(32) 


236(30) 244(21) 


IbJIb 


1.15 


1.10(1) 




1.20(3)(1) 


1.16(4) 1.16 



222.6(16.7) (||) 
266(32) 



Table [H} which have the typical magnitude of the heavy quark corrections. This explains 
the closeness of the obtained values of constants to our previous results [l4J. On the other 
hand, for light mesons, especially for the pion and kaon, the binding energies are not small 
on the meson mass and quark energy scales and, thus, such factor gives no suppression. The 
complete relativistic expression (fTrJj) for decay constants ju brings theoretical predictions 
for light mesons in good agreement with available experimental data. 

The comparison of our values of the decay constants of light mesons with other predictions 
in Table 11111 indicate that they are competitive even with the results of more sophisticated 
approaches (e.g. [llj]) which are based on the Dyson-Schwinger and Bethe-Salpeter equa- 
tions. On the other hand, our model is more selfconsistent than some other approaches 
[IcI llr^ . We calculate the meson wave functions by solving the quasipotential equation in 
contrast to the models based on the relativistic Hamilton dynamics [16( where various ad 
hoc wave function parameterizations are employed. 

In Table HVI we confront our results for pseudoscalar decay constants of the heavy-light 
mesons as well as their ratios with the recent predictions based on the Salpeter equation |l5[ , 
values from the unquenched two- 2] and three-flavour 0, 0| lattice QCD, 6 QCD sum rules 
IE El and available experimental data [l^. Reliable experimental data, up till recently, 
existed only for fo s , which was measured by several experimental collaborations (ALEPH, 
DELPHI, L3, OPAL, Beatrice, CLEO, E653, WA75, BES) both in the D s -> fw and the 
D s — ► tv decay channels. At present, experimental errors are still rather large for this 
constant. Very recently, the CLEO Collaboration published a relatively precise value 
for the decay constant fo measured in D —>■ \w decay. We see from Table ITVl that there is a 
good (within error bars) agreement between all presented theoretical predictions as well as 
with available experimental data. 

In summary, the weak decay constants of pseudoscalar and vector light and heavy-light 
mesons were investigated with the special emphasize on the role of relativistic effects. For our 
calculations we used the meson wave functions which were obtained by the numerical solution 
of the quasipotential equation with the nonperturbative treatment of all spin- dependent 
and spin-independent relativistic contributions to the quark interaction potential. It was 
argued that both positive- and negative-energy parts of the quark propagators in the weak 
annihilation loop should be taken into account. The positive-energy contributions, which are 



6 The recent quenched lattice QCD values |6j] for the pseudoscalar decay constants are fx = 152(6) (10) 
MeV, f D = 235(8) (14) MeV and f Ds = 266(10) (18) MeV. 
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usually considered in the semirelativistic quark models, significantly overestimate the decay 
constants of light mesons. We showed that the negative-energy contributions to the light 
meson decay constants are large and negative. Their account is necessary to bring theoretical 
predictions in agreement with experimental data. On the other hand, these negative-energy 
contributions are considerably smaller for decay constants of heavy-light mesons and have 
the order of magnitude of the lowest correction in the heavy quark expansion. The consistent 
inclusion of relativistic effects coming both from the quark propagators and the meson wave 
functions considerably improve the accuracy and reliability of the obtained predictions. 
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and useful discussions. Two of us (R.N.F. and V.O.G.) were supported in part by the 
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for Basic Research under Grant No. 05-02-16243. 
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